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LIFT FORCE OF AN ARROW-SHAPED WING* 

By M. I. Gurevich 

1. STATEMENT OF THE PROBLEM AND INTRODUCTION 

The flow about a conical body of an ideal compressible fluid, 
is considered. Assume that the velocity of the oncoming flow at 
infinity W is directed, along the z-axis. The system of Car- 
tesian coordinates x, y, z with origin at the vertex of the 
cone 0 is shown in figure 1. From the considerations of the 
dimensional theory, it may be found that along any ray issuing 
from 0 the components of the velocity u, v, W+w along the 
coordinate axes will maintain a constant value. It is further 
assumed that the conical body has such shape and disposition rela- 
tive to the flow that u, v, and w are small in comparison 
with W. The equation of continuity can then, as is known, be 
given in the form 


♦ V - 1) = o (1.1) 

where F may be considered either as the velocity potential cp 
or the components of the velocity u, v, w; M = w/a is the Mach 
number, where a is the velocity of sound at infinity upstream 
of the body. It is assumed that the flow is everywhere super- 
sonic (M >1). 

In view of the linearity of equation (l.l), conical flows 
with different vertices may be superposed on one another. Equa- 
tion (l.l) is of the hyperbolic type, which means that each point 
of the flow may have an effect only on the points located within 
its Mach cone; by a known method, the body may be transformed 
without affecting the flow about Its remaining forward part. It 
must be borne in mind, however, that the Mach cones constructed 
for those points of the body at which the cut is made must 
nowhere Intersect the remaining parts of the body. A. Busemann 
(reference 1) pointed out the analogy that exists between the 
problem of conical flows at supersonic velocity and small addi- 
tional velocities and the two-dimensional problem for subsonic 
flow in the statement of the problem by S. A. Chaplygin (refer- 
ence 2), which becomes exact when the equation of state of the 
gas has the form 


*"0 Pod’emnoi Sile Strelovidnogo Kryla v Sverkhzvukovom 
Potoke." Prikladnaya Matematika i Mekhanika, vol. X, 1946, 
pp. 513-520. 
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where p is the pressure, p the density of the fluid, and 
and C 2 are constants. Making use of the considerations of 

Chaplygin, Busemann reduced the determination of u, v, w within 
the Mach cone to the solution of a Laplace equation in the two 
variables. In sections 2, 3, and 4, a direct derivation is given 
of the formulas of Busemann needed herein and the remaining 
sections are devoted to the investigation of the lift of an arrow- 
shaped wing. 


2. TRANSFORMATION OP EQUATIONS OF CONTINUITY 

For the conical flow, the components u, v, w can be repre- 
sented as functions of the coordinate ratios | ** x/z, tj = y/z. 

The magnitudes | and tj may' be considered as the Cartesian 
coordinates in a plane perpendicular to the z-axis and inter- 
secting it at the distance 1 from the origin (fig. 1). The Mach 
cone for the -po int 0 cuts out in the |,tj plane a circle of > 

radius A = 1 /a/M 2 - 1. 

It is assumed for definiteness that F = w. It is readily 
shown that equation (1.1) in the variables l,r\ will, within the 
circle of radius A, be of the elliptic type and outside it of 
the hyperbolic type. As has already been mentioned, the deter- 
mination of w within the circle reduces to the solution of a 
Laplace equation. This problem outside the Mach cone must be 
separately solved. Both solutions join at the Mach oone in such 
a way that on passing through it the velocity changes continuously. 

Proceed to the transformation of equation (1.1). It is 
reduced to the equation of Laplace in three variables for the 
Independent variables 

x* = ix •JO? - 1 y* = iy - 1 z* = z (2.1) 

In spherical coordinates r, a, and £ where 
x* = r cos a sin •d y* = r sin O sin £ z* = r cos £ 


( 2 . 2 ) 
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equation. (1.1) has the form 

(r 2 F r sin 4) + ^ (jffir?) + |j (F4 sin 4) - 0 

Because of the constancy of the velocity along the rays issu- 
ing from 0, 1 

v a(J + sin A ^ (v^ sin i) *= 0 (2.3) 

Within the Mach cone, imaginary transformat ions may he 
avoided hy setting t - 16; then equations (2.1), (2.2), and 
(2.3) are replaced by 

| = A cos a th 0 i\ a A sin a th 0 (2.4) 

^00+ ® ( W Q S & ®) 13 0 (2.5) 

In order to reduce equation (2.5) to the Laplace equation in 
polar form, there remains to be introduced a new independent 
variable €, determined from the condition de/d0 = e/sh 0; 
whence 


€ = th I v ao + e (w e c) = 0 (2 * 6) 

The same result may be repeated for u and v but not for 
cp because cp r ^ 0. Thus v is the real part of a certain func- 
tion of a complex variable, vhich is denoted by 

Af (e e iCT ) » Af (t) * w + is 

From equations (2.4) and (2.6), it follows that to transform 
the circle with center at 0 and radius A in the plane |,tj 
(interior of the Mach cone) on the unit circle of the t plane 
with center at the origin of coordinates, the radii vectors, keep- 
ing the polar angles a unchanged, must be transformed by the 
formula 


^This condition may be expressed in the form 

xw x + yWy + zw z = 0 whence x^w^ + y*w^ + zi ‘ v z » = rw 2 = 0 
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A th 0 = 


2A th 0/2 
1 + th 2 0/2 


2Ae 
1 + e J 


(2.7) 


t 


3. COMPUTATION CSP VELOCITY COMPONENTS ALONG X- AND T-AXES 

Having tlie function f(T), "by integration the following 
equation may he obtained: 


Cl) » u + iv 


1 

2 



(3.1) 


This equation is obtained as a result of setting up _ 
d(«)=C 0 |d|+a)^ dr] and transforming to the variables t , t. The 

partial derivatives 0)| and co^ are computed with the aid of 

the relations u x = v , u z » w x , v z = w y , and equation (l.l). 

It may be noted Incidentally that according to the analogy of 
Busemann the magnitudes u, v, w correspond to x, y, -\|/ of 
Chaplygin and equation (3.1) agrees essentially with that of 
Chaplygin (reference 2, ch. V, eq. (94)). 


4. BOUNDARY CONDITIONS 


Consider the boundary conditions that are encountered in the 
problem of the arrow-shaped wing. If the wing does not extend 
beyond the Mach cone, then on the circle ( T | * 1 the additional 
velocities u, v, w are equal to 0. 2 Assume a plane wing form- 
ing an infinitely small angle a with the z-axis and an angle 
a with the x-axis. The normal velocity Im(coe" ia ) along It 
will have a constant value. Whence using equation (3.1), the 
following is obtained: 


3m £oe" 




- Ids) 


ds = 0 


+ e(dw + ids) 


-n 


= const 
(4.1) 


2 Along those parts of the Mach cone where dw «* 0, u and v 
will likewise be constant because from equation (3.1) for T = e iCT 

It follows that do) = - i (e ia ids - ids/e -10 ) = 0. 
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Let parts of the arrow-shaped plane wing inclined at an infi- 
nitely small angle to the x, z plane extend beyond the Mach cone 
(fig. 2). The region of influence of the extending parts of the 
wing fills in space the interiors of the Mach cones drawn from each 
point of the wing. In the £ /n plane, this region is represented 
by ACDEEF, which is bounded by the tangents AC, AE, EB, and 
BF drawn from the edges of the wing to the Mach cone and the arcs 
CD and IF. 


The flow outside the Mach cone evidently consists of plane 
flows. Hence, in each of the regions ACG, AGE, BBH, and FHB 
the velocities and the pressures will maintain constant values, 
whereas on passing through the wing w will change sign. For 
a wing symmetrical relative to the axis tj, the boundary condi- 
tions for f(r) on the corresponding arcs of the circle will be 


w = Wq on CG and HD 
w = -Wq on GE and FH 
w = 0 on CD and IF 



(4.2) 


5. ARROW-SHAPED WING EXTENDING BEYOND MACH CONE 

Only a plane arrow-shaped wing with the arrow angle y (fig. l) 
symmetrical with respect to the plane x = 0, forming with the 
plane y = 0 a small angle 8 and finally cut along a straight line, 
which shall be assumed approximately coinciding with the axis n = 0, 
is considered. 


The analysis begins with the case where the ends of the wing 
extend outside the Mach cone (tan y > A). The function f(T) can 
easily be oonstruoted from its singularities. 


w + is = Af ( T) 


i(0 0 

— ** 


T 2 -e~ 2iq b 

T 2 e- 2i<7 0-l 


(5.1) 


It is not difficult to be convinced of the validity of equa- 
tion (5.1) by direct check: of conditions (4.1) and (4.2) because 
the analytic character of Af(T) within the circle |t| = 1 is 
evident . 


The concepts w Q -.and Cq must be expressed in terms of the 

given magnitudes, the angle of attack 8 (or, what amounts to the 
same thing, in terms of the vertical velocity on the wing Vq = -W8), 

and the cone angle y. From figures 1 and 2 it is readily seen that 
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tan y = 


008 o 


0 


(5.2) 


The velocity v Q may he expressed in terms of w Q hy inte- 
gration of equation (3.1): 


v o - ( T - ^)f^ dT 


v 0 sln g c 

JtA 


Im 


Ig -(T-e- lg O)( T -e lg O ) 
(T+e" ia O)( x+e iCT 0) 


“ 10 ' 


whence 


AVr 


Vr 


sin O n 


or 


^ = flA 
y sinCT. 


(5.3) 


The lift-force coefficient of the wing is now computed. 



2Y 

PW 2 S 



7 

w d| 
7 


(5.4) 


where Y is the lift force, p the density of the gas, and 
S = tan 7 is the area of the wing. Noting that along the parts of 
the wing extending outside the Mach angle oone, w = ±w Q , the fol- 
lowing is obtained: 


C Y 


2 

W tan y 


2wg(tan y - A) 



or integrating hy parts 


i 

i 

c 

■for] _ A 

r 1 t 

d(w+is) 

J ~ W tan y 

7 1 

! l+T 

• w 0 


J 

-1 



(5.5) 


/ 
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and finally 


4w 0 

4A 

tan y - sin 20 

S 

f— 1 

W tan y 


-1 1 + T 2 (T 2 -e- 2ia 0 )(T 2. e 2ia 0) 1 


(5.6) 


The integral entering equation (5.6) is most conveniently 
computed by deforming the contour of integration into an upper semi- 
circle of unit radius passing around the poles = e-*-°b, 

T 2 = and = i over infinitesimally small semicircles. 

From equation (5.5) and the boundary conditions on the Mach cone, it 
is easy to see that the integrals over the arcs of the upper semi- 
circle |t| = 1 will be imaginary and evidently cancel each other. 
The computation of thus reduces to the finding of the half- 

residues at the points T^, Tg, and T^. By computation 


pl 


T a l rafcfll * _2 — : . tan o r 


l-i 


l+T* 


V r 


COS C„ 


whence making use of equations (5.6) and (5.3), 


C T « 4 PA = =H - (5.7) 

This equation agrees with the well-known formula of Ackeret 
(reference 3) for the plane wing, which is a particular case of the 
arrow-shaped wing (for y = k/2). 

The distribution of the lift force, or more accurately of the 
magnitude w/(AWP) proportional to the intensity of the lift force 
for the arrow-shaped wing extending outside the Maoh cone, is shown 
in figure 3 (curve II). On the parts of the wing outside the Mach 
cone, the intensity of the lift force is constant. Within the Mach 
cone, curve II was computed by equation (5.8) obtained from eq.ua- 
tion (5.1): 
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v 

AWP 


sin O r 


1 - 1 arc tan 22L?P IzL l 


jt 


sin <3, 


'O 1+T 


(5.8) 


On the axis of ahscissas the values of i referred to ha lf 
the wing span were laid off. The Mach cones in figure 3 are rep- 
resented in the form of the dotted semicircles. The nearer the 
Mach cone approaches the edges of the wing, the more intense is the 
lift force at the edges (for tan y/k = 1, curve III, at the edges 
of the wing w/ (A WP) =»). For A/tan y = 0 (case of Ackeret), 
the intensity of the lift force along the wing is constant. Curve IV 
gives the distribution of the lift-force intensity for a wing within 
the Mach cone. (See section 6.) In this case at the edges of the 
wing, w/(AWP) = 00 . 


6. ABHOW-SHAPED WING LOCATED WITHIN MACH CONE 

The case where the wing does not extend beyond the Mach cone 
(tan y < A) is now considered. In the plane T, the wing is rep- 
resented by the segment of the real axis from -b to b where 
according to equation (2.7) 


tan y = — (6.1) 
1+b* 2 

Conditions (4.1) and (4.2) in the plane t assume the form: 
w = 0 on the circle, and s = 0 on the wing. 

From the symmetry of the boundary conditions for Af (T), it 
follows* that in passing through the wing w changes sign. The 
function Af(t) may be continued by the principle of mirror reflec- 
tion over the entire plane with cuts from -oo to -l/b, from 
-b to b, and from l/b to « (fig. 4). It is simplest to con- 
struct w+is = Af(T) from the singular points and the zeros: 


w+is = B 


T 2 +l 


Aj^-^Hi/b^T 2 ) 


Af 


( 6 . 2 ) 


where B denotes a real constant to be determined. Equation (6.2) 
can also be checked directly. The value of Vq is found on the 

plate. For this value, the integration in equation (3.1) is carried 
out over the imaginary radius. Because along the axis of imaginaries 
T = -T and -df/dT = df/dT; hence 
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r 0 


^T=0 = " 2 


ifzi « dT - - B O-tl/b)* 
T dT 2A 


' a -/) 8 d T 

0 (i 2 -i 2 ) 3/2 dA 2 -T 2 ) 3/2 


whence 


_ B(b+l/b)‘ 


^1 


2A 


(1+e 2 ) 2 de 


0 (b 2 +e 2 )(l/b 2 +e 2 ) Nj(b 2 +e 2 )(l/b 2 +e 2 ) 


(6.3) 


setting 

€ « b tan <P / Jl-b^ » k b 2 » k' tan cp^ = l/ /^k’ 

(6.4) 

the integral on the right aide of equation (6.3) is readily trans- 
formed to the canonical form. 

^ B.-(l-k») sin 2 <p] 2 dtp 

_ Q (1-k 2 sin 2 cp) aJ 1-k 2 sin 2 qp 

or 


T 0 " 


_ Bb (b+l/b)‘ 


2A 


v 0 = - 


B 


2A A/k* 


f*l 


d«p 


(1-k 2 sin 2 cp) /^1-k 2 sin 2 <p 


+ Ek^C^lO+BC^lc) 


where F(*¥> 1 ,k) and B (<?■]_, k) are elliptic integrals of the first 
'and second kind. As is known (ref. 4, equation (126)) 


rfPl' 

dq> 

J ° (1-k 2 e in2« 3/2 


B«P x ,k) 

k’ 2 


fc 2 Bin <P-j_ cos qp^ E(Cp 1# k) i-k» 

k l2 /\J 1-k 2 sin 2 ^ k‘ 2 k ' 2 


whence 

B = 4k 7 

At 0 = kT^k) + Kcp^k) - (l-k*)/2 • 


(6.5) 
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With the aid of equation (6.5) from equation (6.2), and from 
equation (5.1) with the aid of equation (5.3), it is found that in 
the limiting case when the •sing touches the Mach cone (o = 0, 


b - 1) 


w+is 




T c +1 

T 2 -! 


( 6 . 2 a) 


The lift-force coefficient of the wing is now computed. From 
equation (5.4) transforming with the aid of equation (2.7) to the 
variable € , 


.. 

4 A 

W tan y v 

J-t 


l-i 2 

(l+e 2 ) 2 


dc 


( 6 . 6 ) 


The integral entering equation ( 6 . 6 ) may be replaced by a 
contour integral about the wing. 


r \ 

2A 

W tan 7 c 


(w+ie) 


1-T 2 

(1+T 2 ) 2 


dT 


(6.7) 


This integral is most simply computed by deforming the contour 
of integration into a unit circle passing around the poles at the 
points T = dki over infinitesimally small semicircles. Because 
on the cirole w is equal to 0 , it is easily seen that the inte- 
grals over the arcs of the circle will be imaginary and will mutually 
cancel. The computation of Cy thus reduces to finding the semi- 

residues at the points T = ±i. 


2AB 
W tan 


1-T 2 dT = 2*B 

/^(b 2 -T 2 )(l/b 2 -T 2 ) 1+T 2 W 


whence 


£z = 2 it *Iis? 

PA + E^k) - (l-k' )/2 


( 6 . 6 ) 
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From e quations (5.7), (6.1), and (6.8), it is seen that 

CyAjM 2 - 1/(43) is a function of aJm 2 - 1 tan 7. The computations 

(fig. 5) can he carried out by assigning various values to the 
parameter 3c. The center of pressure of the wing is at the center 
of gravity of the triangle. 


7. RMABKS ON DRAG 


If the wing extends beyond the Mach cone, the velocities at 
the edges of the wing will be finite and the efficiency of the 
wing Y/x * l/p, where X is the drag. In the case where the 
wing Is entirely within the Mach cone, then as a result of the 
infinite velocities at the leading sharp edges suction forces may 
be expected, which increase the efficiency of the wing. 


The drag of the wing can be computed with the aid of the momen- 
tum theorem applied to the volume of gas enclosed in the cylinder 
about the wing and the Mach cone, one of the bases of the cylinder 
lying in the £,t] plane and the other to the left of the origin of 
coordinates in the region of the undisturbed flow.. 


X 


P 

2 \j 


[u 2 + v 2 + (M 2 -l) w 2 )^] d| dti 


(7.1) 


where the integration Is extended over the entire base of the 
cylinder in the | ,tj plane. In deriving equation (7.1), use is 
also made of the Bernoulli integral, in which after expanding in 
powers of u, v, w all terms were neglected that contained the 
additional velocities to degrees higher than the second, the adia- 
batic condition, and the equation of continuity in integral form. 


Translated by S. Reiss, 
National Advisory Committee 
for Aeronautics 
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